OVERTONES OF ISOSCALAR GIANT RESONANCES IN MEDIUM-HEAVY 
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A semi-microscopic approach based on both the continum-random-phase-approximation (CRPA) 
method and a phenomenological treatment of the spreading effect is extended and apphed to describe 
the main properties (particle-hole strength distribution, energy-dependent transition density, partial 
direct-nucleon-decay branching ratios) of the isoscalar giant dipole, second monopole, and second 
quadrupole resonances. Abilities of the approach are checked by description of gross properties 
of the main-tone resonances. Calculation results obtained for the resonances in a few singly- and 
doubly-closed-shell nuclei are compared with available experimental data. 
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I. INTRODUCTION 



Experimental and theoretical studies of high-energy gi- 
ant resonances (GRs) have been undertaken in recent 
years to understand better how the different characteris- 
tics associated with GR formation (concentration of the 
particle-hole strength, coupling to the continuum, the 
spreading effect) are affected with increasing the GR en- 
ergy. Many of known high-energy GRs are the next vi- 
bration modes (the overtones) relative to the respective 
low-energy GRs (the main tones). The lowest energy 
overtone is the isoscalar giant dipole resonance (ISGDR) 
experimentally studied in a few medium-heavy and heavy 
nuclei via the {a, a') reaction 0,111111. The ISGDR is 
the overtone of the 1~ zero-energy spurious state (SS), 
associated with center-of-mass motion. The experimen- 
tal results of Refs. Q, H H Q are concerned with distri- 
bution of the respective dipole strength. Only recently 
direct nucleon decays of the ISGDR have been firstly 
observed using the {a,a'N) reactions M H 01- These 
studies are planned to be continued Q- The isovec- 
tor giant charge- exchange (in the /3~ channel) monopole 
and spin-monopole resonances are the overtones of the 
isobaric analogue and Gamow- Teller resonances, respec- 
tively. These overtones have been studied via charge- 
exchange reactions H ^3 ■ Direct proton decays of the 
giant spin-monopole resonance have been recently ob- 
served using the (^}ie,tp) reaction ^^l- Other candidates 
for studies of high-energy GRs are the overtones of the 
isoscalar giant monopole and quadrupole resonances (IS- 
GMR2 and ISGQR2, respectively). The corresponding 
main tones, having relatively low energy, have been much 
experimentally studied 0,|1I1]. However, the evidence 
for existence of the ISGQR2 has apparently been found 
only recently H |1| . 

Theoretical studies of the isoscalar overtones deal pri- 
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marily with the ISGDR. Microscopically, this GR is 
mainly due to SHlu particle-hole- type excitations, while 
the overtones of the isoscalar quadrupole and monopole 
GRs are mainly due to Ahuj excitations. First micro- 
scopic calculations of the strength distribution for the 
ISGDR and for all the above-mentioned isoscalar over- 
tones have apparently been done in the eighties (Refs. pl| 
and (T^, respectively). In recent years Hartree-Fock -|- 
RPA calculations with the use of the Skyrme interactions 
have been done to specify the ISGDR energy and, after 
comparing with experimental data, to make a conclu- 
sion about the nuclear incompressibility [isl IT^ . Similar 
goals are pursued in the approaches based on the rel- 
ativistic version of the RPA and on the semi-classical 
treatment of nuclear vibrations (see, e.g., Refs. and 
|16|. respectively). However, the microscopic structure 
and "differential" properties of high-energy GRs, corre- 
sponding to collective excitations of nuclei as finite-size 
open Fermi-systems, are of particular interest. Attempts 
to describe the main properties (the strength distribu- 
tion, energy-dependent transition density, partial direct- 
nucleon-decay branching ratios) of the ISGDR have been 
undertaken in Refs. [TtLIi^ within a semi-microscopic ap- 
proach, based on both the CRPA method and a phenome- 
logical treatment of the spreading effect. In these refer- 
ences the gross properties (parameters of the strength 
distribution, transition density) have been satisfactorily 
described for the ISGDR in a few medium-heavy and 
heavy nuclei. Contrary to the gross properties, the par- 
tial branching ratios for direct nucleon decay of a partic- 
ular GR carry information on its microscopic structure 
and also on coupling to the continuum and the spread- 
ing effect. The attempts (undertaken in Refs. 0, 0| 
for the ISGDR in ^°*Pb) to evaluate the partial direct- 
nucleon decay branching ratios met with difficulties con- 
cerned with taking the spreading effect into account. The 
way to overcome these difficulties has been drawn in 
Ref. The overtone of the isoscalar monopole res- 

onance was theoretically studied within a CRPA-based 
approach in Ref. poj mainly for searching for narrow 
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("trapped") resonances, having a small relative strength. 
The gross properties of the ISGMR2 are described in 
Ref. !21|, where the microscopic Hartree-Fock + RPA 
approach with the use of the Skyrme interactions and 
also a semi-classical approach are employed. Within the 
semi-microscopic approach, the main pro perties of this 
resonance are briefly described in Ref. [l9|. 

Motivated by both aspirations to describe the main 
properties of high-energy GRs and forthcoming experi- 
mental data concerned with the isoscalar overtones, we 
pursue the following goals in the present work: 

(i) extension of the CRPA-based semi-microscopic ap- 
proach to describe direct nucleon decay of high- 
energy GRs; 

(ii) description of the gross properties of the isoscalar 
main-tone resonances to check abilities of the ap- 
proach and to find out probing operators appropri- 
ate for overtone studies; 

(iii) calculation of the partial dircct-nucleon-decay 
branching ratios for the high-energy component of 
the ISGDR; 

(iv) description of the main properties of the sec- 
ond isoscalar giant nionopole and quadrupole reso- 
nances; and 

(v) comparison of the calculation results, obtained for 
singly- and doubly-closed-shcU nuclei ^^Ni, ^"Zr, 
"'^^^Sn, ^''^Sm, and ^"*Pb, with available experimen- 
tal data. 

The paper is organized as follows. In Sect. |n] the ba- 
sic elements, ingredients and new points of the approach 
are presented. Sect. Illll contains calculation results and 
available experimental data on properties of the isoscalar 
giant main-tone and overtone resonances. Discussion of 
the results and concluding remarks are given in Sect. lIVI 



II. BASIC ELEMENTS AND INGREDIENTS OF 
THE APPROACH 

The continuum-RPA method and a phenomenologi- 
cal treatment of the spreading effect are the basic ele- 
ments of the semi-microscopic approach. In implemen- 
tations of the approach the following phenomenological 
input quantities are used: a realistic nuclear mean field 
and the Landau-Migdal partical-hole interaction bound 
together by selfconcistency conditions; an energy- and 
radial-dependent smearing parameter. 



A. CRPA equations 

The CRPA equations are taken in the form accepted 
within the Migdal's finite Fermi-system theory As 



applied to description of isoscalar particle-hole-type ex- 
citations in closed-shell nuclei, these equations are given 
in detail in Refs. fvA ITsI l and not shown here. (For 
reader's convenience we use (in the main) the notations 
of Refs. (17nl8i] and sometimes refer to equations, tables 
and figures from these references). The nuclear polariz- 
ibility Pl(w), strength function S'l(w), and the energy- 
dependent transition density pL(r, w) = pL{r,u>)YLMin) 
(lu is the excitation energy) , corresponding to an isoscalar 
probing operator Vt(r) = VL{r)YLM{n), can be cal- 
culated to describe gross properties of the respective 
isoscalar GR within the CRPA. The above-listed quanti- 
ties are expressed via the radial part of the effective prob- 
ing operators V^(r, w) (a — n,p is the isotopic index), 
which are different from VL{r) due to core-polarization 
effects caused by the particle-hole interaction. The ui- 
dependence of the above quantities comes from that of 
the free particle-hole propagator. The latter can be 
expressed in terms of occupation numbers n'^^, radial 
bound-state single-particle wave functions and ra- 
dial Green functions g"x)i^^ ^' ^ ^ij.^'^) to take exactly the 
single-particle continuum into account (see Eqs. (l)-(3) 
of Ref. [13 and Eq. (1) of Ref. [3). 

When compared with Ref. ^2^, the new element of 
the CRPA equations is used within the approach — the 
direct-nucleon-escape amplitude M^(tt') (see Eqs. (4), (5) 
of Ref. 0|]). This amplitude is proportional to the prod- 
uct of {n'jiy^^ and the matrix element of the respective 
radial effective operator taken with the use of the radial 
bound-state wave function and the radial continuum- 
state wave function x[^l" (m is the set of quantum num- 
bers for an occupied single-particle level; e = + a; and 
(A) are the energy and quantum numbers of an escaped 
nucleon, respectively, c = fi, (A), a is the set of nucleon- 
decay-channel quantum numbers compatible with the re- 
spective selection rules) . The w-dependencc of the direct- 
nucleon-escape amplitude comes not only from the effec- 
tive probing operator but also from the continuum-state 
wave function. 

The partial direct-nucleon-decay branching ratio b^{d) 
can be reasonably defined as the ratio of the squared 
nucleon-escape amplitude integrated over a certain 
excitation-energy interval (5 = wi — 072 to the respective 
strength function integrated over the same interval. The 
total branching ratio fo^j = J^c is equal to unity for ar- 
bitrary interval i5, as it follows from the unitary condition, 
which is valid within the CRPA (Eq. (4) of Ref. [13). 

For comparison with the experimental branching ratios 
we sometimes replace the occupation numbers in the 
expression for 6^ = with the respective experi- 

mental spectroscopic factors 5'^. In such a way we take 
phenomenologically into account coupling of single-hole 
states fi~^ populated after direct nucleon decay of GRs 
to low-energy collective states. 

In description of the gross properties of high-energy 
GRs, nucleon pairing in open-shell subsystems can be 
neglected with a high accuracy. To take into account 
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TABLE I: The peak energy, total width (both in MeV), and 
parameters rjL (in fm^) calculated for the L = 0, 1, 2 isoscalar 
GRs in nuclei under consideration. The respective experimen- 
tal values (riven with errors) are taken from Ref. 0] (^ = 0, 1) 
and Ref. ^ {L = 2). 





s^Ni 


*Zr 


"«Sn 




20«Pb 


ISGMR 














17.7 


16.2 


15.7 


14.8 


13.6 






16.9 ±0.1 


15.4 ±0.1 


15.4 ±0.1 


13.5 ±0.2 


r 


4.4 


4.2 


4.0 


3.9 


3.9 






4.3 ±0.1 


5.7 ±0.3 


3.9 ±0.2 


4.0 ±0.4 


ISGQR 














14.7 


12.7 


11.9 


11.2 


10.0 






14.0 ±0.2 


13.2 ± 0.2 


12.2 ±0.2 


11.0 ±0.2 


r 


2.8 


2.65 


2.6 


2.4 


2.2 






3.4 ±0.2 


3.3 ±0.2 


2.4 ±0.2 


2.7 ±0.3 


ISGDR 












rjL=i 


22.71 


29.28 


34.74 


39.74 


51.07 


LE 
peak 


9.2" 


12.1 


9.3 


12.0 


7.0' 






16.8 ±0.4 


15.5 ± 0.4 


13.0 ±0.3 


13.0 ±0.2 


^peak 


27.4 


26.0 


26.1 


25.0 


22.9 






26.9 ±0.7 


24.9 ± 0.7 


25.0 ±0.3 


22.8 ± 0.3 


ISGMR2 












f?L=0 


42.04 


52.75 


55.57 


54.50 


75.25 


HE 
^peak 


30.8 


29.9 


32.0 


33.8 


32.1 


ISGQR2 












r]L=2 


24.38 


30.50 


35.90 


40.41 


52.54 


^.HE 

^peak 


32.9 


34.6 


34.1 


32.8 


30.5 



"The LE-component has additional maximum at 12.1 MeV 
''The LE-component has additional maximum at 11.4 MeV 



the effect of nucleon pairing on the direct-nucleon-decay 
branching ratios in the CRPA equations we replace the 
occupation numbers with the respective Bogolyubov's 
factors ti^. The latter can be calculated in an isospin- 
selfconsistent way with the use of experimental pairing 
energies |23|. To take phenomenologically into account 
coupling of singlc-quasi-particle states populated after di- 
rect nucleon decay to low-energy collective states the cal- 
culated 6^ value is divided by i;^ and multiplied by S*^ 

M 

An important aspect of theoretital studies of GR 
overtones within the RPA is the choice of an appro- 
priate probing operator. It is convenient to choose 
this operator with the condition that the main tone is 
not being excited. In this case the overtone exhausts 
most of the respective particle-hole strength. As ap- 
plied to description of the ISGDR, charge-exchange gi- 
ant monopole, and spin-monopole resonances the choice 
of appropriate probing operators is discussed in Refs. 0, 
l2^ . and [2j|, respectively. In particular, the radial 
part of the isoscalar second-order dipole probing op- 
erator V^^\(r) = — rjL^ir is used for description 



TABLE II: Paramaters of the ISGMR and ISGQR calculated 
for a certain excitation-energy interval. All the parameters 
are given in MeV except for x, which is given in %. 



rmcl. 


Ull — UJ2 


ISGMR 


A 


X 


UJl — UJ2 


ISGQR 


A 


X 




12-31 


18.21 


3.06 


89 


10-20 


14.89 


1.90 


70 


9»Zr 


10-25 


16.47 


2.64 


87 


9-18 


13.07 


1.78 


70 


"«Sn 


10-20 


15.38 


2.13 


79 


8-19 


12.51 


2.13 


72 


i^^Sm 


10-20 


14.85 


2.05 


80 


7-19 


11.91 


2.31 


74 


2°«Pb 


10-20 


13.89 


2.05 


79 


6-17 


10.51 


2.18 


71 



of the ISGDR. Parameter ?7l=i is defined by the con- 
dition: / Pl=i(^)^l=i('^)''^'^^ — 0, where /o|£i(?') is 
the spurious-state transition density. To describe the 
properties of the second isoscalar giant monopole and 
quadrupole resonances the radial part of the respective 
second-order probing operators is taken in the form: 
vj^^\r) = — r]Lr^. The parameters 77^ in this expres- 
sion are determined by the condition: 

j PL{r,ojpeak)VP{rydr = Q, (1) 

where PhifT^peak) is the radial part of the energy- 
dependent main-tone transition density taken at the peak 
energy of the main-tone strength function. 

B. Smearing procedure 

Within the approach, the speading effect, i.e. cou- 
pling of particle-hole-type doorway states to many-quasi- 
particle configurations, is phenomenologically taken into 
account in terms of a proper smearing parameter. Some- 
what different smearing procedures are used for descrip- 
tion of low- and high-energy GRs. The above-mentioned 
cij-dcpcndent quantities calculated within the CRPA for 
a low-energy ( "subbarrier" ) GR can be expanded in 
terms of the Breit-Wigner-type doorway-state resonances 
nonoverlapped on their total escape widths F^. To take 
the spreading effect in evaluation of the energy-averaged 
(jj-dependent quantities into account, the doorway-state 
resonances are smeared independently of one another. 
The smearing procedure comprises the replacement of 
by PT + /, or (what is the same) the replacement of 
w by a; ± i//2. The parameter /, the mean doorway- 
state spreading width, is fitted to reproduce the experi- 
mental total width of the considered GR in calculations 
of the energy-averaged strength function. To calculate 
the energy-averaged direct-nucleon-escape amplitudes it 
is also necessary to average the potential-barrier penetra- 
bility in a proper way because the directly escaped nu- 
cleons have relatively low ("subbarrier") energies. Thus, 
the energy-averaged transition density and partial direct- 
nucleon-decay branching ratios can be evaluated within 
the approach without the use of any free parameters. 
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TABLE III: Parameters of the ISGDR calculated for different 
excitation-energy intervals. All the parameters are given in 
MeV except for x, which is given in %. 



LE-ISGDR 



HE-ISGDR 



nucl. 


LOl — UJ2 




A 


X 




LO 


A 


X 




5-16 


11.56 


2.08 


12 


16-36 


26.62 


4.55 


75 












16-60 


27.94 


6.36 


86 


9»Zr 


11-18 


13.78 


2.10 


10 


18-32 


25.42 


3.48 


68 




5-16 


11.32 


2.19 


11 


16-40 


25.94 


4.74 


81 












16-60 


26.69 


6.14 


86 


i«Sn 


11-18 


13.91 


1.99 


11 


18-32 


25.09 


3.23 


67 




5-15 


10.50 


2.36 


13 


15-35 


24.78 


4.22 


75 












15-60 


25.91 


6.13 


84 


"*Sm 


5-15 


10.74 


2.22 


12 


15-35 


24.20 


4.04 


77 












15-60 


25.22 


5.92 


85 


208p^, 


8-15 


11.18 


1.89 


14 


15-24 


20.73 


2.39 


42 




5-15 


9.90 


2.55 


17 


15-35 


23.02 


3.90 


74 












15-60 


23.96 


5.83 


81 



This method has been employed in Ref. to decribe 
quantitatively direct nucleon decay of the ISGMR in a 
few nuclei. As applied to direct nucleon decay of the 
ISGDR in ^°^Pb, the method results only in qualita- 
tive description because the respective doorway-state res- 
onances are significantly overlapped 0|. The above- 
outhned smearing procedure is based on a statistical 
assumption: after energy-averaging the doorway states 
"decay" into many-quasi-particle configurations indepen- 
dently of one another . Such an assumption seems to 
be reasonable in view of the complexity of many-quasi- 
particle configurations at high excitation energies. 

As follows from CRPA calculations of the high-energy 
GR strength function, the doorway-state resonances are 
overlapped. In such a case, the smearing procedure 
comprises the replacement of w by w -I- il /2 in the 
CRPA equations. This replacement implies, in fact, 
the use of the imaginary part of the single-particle 
potential, =F(i/2)/(r, w), when the radial Green func- 
tions and continuum-state wave functions are calcu- 
lated. As a result, the energy-averaged w-dependent 
quantities (strength function Sl(!jj), transition density 
PL{r,uj), direct-nucleon-escape amplitude M^iuo)) can 
be calculated at once with the use of a radial- and 
energy-dependent smearing parameter in the CRPA 
equations. In accordance with the statistical assumption, 
the following parameterization for / is used: /(r, oj) = 
I{'^)fws{r, R* ,0.), where f^s is the Woods-Saxon function 
taken with R* > R (R and a are the radius and diffuse- 
ness of the isoscalar part of the nuclear mean field, respec- 
tively). The calculated strength distribution and transi- 
tion density of high-energy GRs are found to be almost 
independent of the "cut off" radius at R* > (1.7 — 1.9)i?. 
The energy-dependent part of the smearing parameter 
/(w) is taken in the form: 




15 20 

0) (MeV) 

FIG. la: The calculated relative energy-weighted strength 
function for the ISGMR. The thick, thin, dashed, dotted, and 
dash-dotted lines are for ^^m, ^°Zr, "®Sn, ""'Sm, and 2°«Pb, 
respectively. 
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FIG. lb: The reduced energy-dependent transition density 
calculated at the peak energy of the ISGMR (Table HJ. The 
notations are the same as in FIG. la. 



I{lu) = a 



(c.-A)2 



1 + (CJ- A)2/B2' 

0, w < A 



w > A; 



(2) 



with universal parameters. Such an energy dependence 
is used for the absorption potential in some versions 
of the optical model for nucleon-nucleus scattering [25l |. 
The above-outlined smearing procedure has been used 
in Ref. 18] to describe quantitatively the gross proper- 
ties of the ISGMR and ISGDR in a few nuclei. How- 
ever, this procedure can not be directly applied to eval- 
uation of the energy-averaged escape amplitudes M^(w) 
(and, therefore, the respective branching ratios) because 
of nonphysical absorption of escaped nucleons outside 
the nucleus. For this reason, in the present work all 
the energy-averaged w-dependent quantities are calcu- 
lated using the "cut off" radius R* = R together with 
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TABLE IV: Comparison of calculated and experimental parameters of the ISGMR and ISGDR. The experimental data (given 
with errors) are taken from Refs. |2fll | and 0, respectively. All the parameters are given in MeV. 







''"Zr 




"«Sn 




208pj3 




ISGMR 




17.89 ± 0.20 


16.5 


16.07 ±0.12 


15.4 


14.17 ±0.28 


13.9 




A 


3.14 ±0.09 


2.6 


2.16 ±0.08 


2.1 


1.93 ±0.15 


2.1 


LE-ISGDR 


Lb 


16.2 ± 0.8 


13.8 


14.7 ±0.5 


13.9 


12.2 ±0.6 


11.2 




A 


1.9 ±0.7 


2.1 


1.6 ±0.5 


2.0 


1.9 ±0.5 


1.9 


HE-ISGDR 


Lb 


25.7 ±0.7 


25.4 


23.0 ±0.6 


25.1 


19.9 ±0.8 


20.7 




A 


3.5 ±0.6 


3.5 


3.7 ±0.5 


3.2 


2.5 ±0.6 


2.4 



0.1 



0.0 




10 12 14 16 

m (MeV) 




' ' I ' ' ' 1 1 ' ' ' 1 1 ' ' 1 1 1 ' ' I ' I ' 1 1 ' I ' 1 1 ' I ' I ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

1 23456789 10 11 12 



r(im) 



FIG. 2a: The calculated relative energy-weighted strength 
function for the ISGQR. The notations are the same as in 
FIG. la. 



FIG. 2b: The reduced energy-dependent transition density 
calculated at the peak energy of the ISGQR (Table 0. The 
notations are the same as in FIG. la. 



the properly increased value of the intensity a in Eq. 1^)1. 
Such a choice allows us to describe correctly the respec- 
tive single-particle resonances in the energy dependence 
of the continuum-state wave function, while the calcu- 
lated gross properties of a high-energy GR are found 
almost the same within both versions of the smearing 
procedure. In Refs. 0, ]|2^ the branching ratios cal- 
culated for a few overtones in ^°®Pb are overestimated 
because of ignoring absorption of escaped nucleons inside 
the nucleus. This shortcoming is partially eliminated in 
Ref . J^J , where all w-dependent quantities are calculated 
using the smearing parameter with intermediate radius 
R* = 1.3i? to take the spreading effect into account in 
evaluation of the branching ratios for the ISGDR in a few 
nuclei. 

In conclusion of this Subsect., we define the energy- 
averaged quantities suitable for description of the 
GR main properties. The relative energy- weighted 
strength function yL{^) = ujSl{(^)/ {EWSR)l is used 
to show exhaustion of the respective energy-weighted 
sum rule (EWSR) by a particular GR. The use 
of the reduced energy-dependent transition density 
i?L(r, w) — r'^ pL{r,u!)S^^^^ (uj) normalized by the con- 
dition / RL{r,u!)VL{r)dr = 1 is convenient to compare 
the transition densities related to different energy re- 



gions jl8l | . The squared and properly normalized energy- 
averaged direct-nucleon-escape amplitude determines the 
respective differential partial branching ratio: 



dbjlicb) 



cLu 



SLioj)dui 



(3) 



These quantities show how the partial and total branch- 
ing ratios are formed. 



C. Ingredients of the approach 

Within the approach, a phcnomcnological isoscalar 
part of the nuclear mean field (including the spin-orbit 
term) and the (momentum-independent) Landau-Migdal 
particle-hole interaction are used as the input quantities 
for CRPA calculations. Parameterization of these quan- 
tities is explicitly given in Ref. Within the RPA, 
the isospin symmetry of the model Hamiltonian can be 
restored. As a result, the isovector part of the mean 
field is calculated selfconsistently via the Landau-Migdal 
isovector parameter /' and the neutron-excess density 
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m (MeV) 

FIG. 3a: The calculated energy-averaged strength function 
for the ISGDR. The notations are the same as in FIG. la. 



TABLE V: Comparison of calculated and experimental pa- 
rameters of the ISGMR and ISGQR in '^^Ni. The experimen- 
tal data are taken from Ref. fsof . All the parameters are given 
in MeV. 



ISGMR 




20.30lJ;?^ 


18.2 




A 




3.1 


ISGQR 


Uj 


16.1 ±0.3 


14.9 




A 


2.4 ±0.2 


1.9 



(see, e.g., Ref. [25). The mean Coulomb field is calcu- 
lated also selfconsistently via the proton density. 

Within the CRPA, the distribution of the isoscalar 
dipole strength (corresponding to the probing operator 
with the radial part VL=\(r) = r) can be calculated for 
a given model Hamiltonian. In studies of Refs. [T^ ITsj. 
the Landau-Migdal isoscalar parameter is chosen for 
each nucleus to make the energy of the 1~ spurious state 
close to zero. If the translation invariance of the model 
Hamiltonian was fully restored, the spurious state would 
exhaust 100% of the EWSR corresponding to the above- 
mentioned probing operator. Within the current version 
of the approach, the SS exhausts in all cases more than 
92% of this sum rule, the exact percentage being depen- 
dent on the nucleus [13, [ill • The small part of the spuri- 
ous strength (less than 8%) is distributed mainly among 
isoscalar dipole Ihuj particle-hole-type excitations. 

The calculation results presented in Sect. IIIll are ob- 
tained with the use of model parameters taken, in the 
main, from previous studies of Refs. fvA IT^ . The mean- 
field parameters and the Landau-Migdal isovector param- 
eter /' = 1.0 are taken from Ref. where the experi- 
mental nucleon separation energies have been satisfacto- 
rily described for closed-shell subsystems in nuclei with 
A — 90 — 208. The isoscalar Landau-Migdal paramater 
/™ is taken equal to 0.0875 in agreement with the sys- 
tematics of Refs. 11^,!!^, while the values of parameter 
/'^^ = -(2.7 - 2.9) are found in Ref. [H in the above- 




to (MeV) 



FIG. 3b: The calculated relative energy- weighted strength 
function for the ISGDR in ^^m. 



described way for each nucleus under consideration. The 
universal parameters of Eq. (01 used for description of the 
spreading effect a = 0.125 MeV'^, A = 3 MeV, B = 7 
MeV are taken the same, as in Ref. e xcep t for the 
a value (a = 0.085 MeV"^ is used in Ref. together 
with R* = 1.8R). 



III. PROPERTIES OF THE ISOSCALAR GIANT 
MULTIPOLE RESONANCES 

A. Gross properties of the main-tone resonances 

We start with description of the main-tone resonances. 
This description allows us 

• to check the quality of restoration of the translation 
invariance of the model within the CRPA; 

• to check abilities of the semi-microscopic approach 
because the gross properties of the isoscalar gi- 
ant monopole and quadrupole resonances are ex- 
tensively studied experimentally; 

• to evaluate the parameters rjL in the expression for 
the radial part of the second-order probing opera- 
tors appropriate for microscopic studies of the over- 
tones. 

The 1" spurious state associated with the center-of- 
mass motion is the lowest energy main-tone state. The 
method to find out properties of the SS within the CRPA 
is described in detail in Ref. 0| (see also Subsect. Ill C|l . 
The method allows us to specify the value of the Landau- 
Migdal parameter f^^ and to calculate the character- 
istics of the SS, relative isoscalar dipole strength xf£i 
and the transition density pf£i(?'). The calculated val- 
ues of and a;f£]^ for nuclei under consideration (ex- 
cept for 58Ni) are given in Table I of Ref. [H. For ^^Ni 
we obtain f'' = -2.646 and xfi^ = 96%. The values 
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FIG. 3c: The reduced energy-dependent transition density 
calculated at the peak energy of the HE-ISGDR (Table 0. 
The notations are the same as in FIG. la. 

= -2.789 and xfi-^ = 94% are obtained for "^Sn 
with taking neutron pairing into account. The radial 
dependence of the calculated transition density pf£i(r) 
is found to be close to dp{r)/dr, where the ground-state 
density p{r) is determined by the bound-state wave func- 
tions for the occupied levels: 

E j:^^^K(x-ir)r. (4) 

a—n,p ^ 

This result is quite expectable, because the calculated 
relative strength a;|£j^ is rather close to 100%. For the 
same reason, the parameter rjL=i in the expression for 
the second-order dipole operator (Subsect. lII Ap is almost 
equal to the value 5 < > /3 (averaging is performed 
over the ground-state density of Eq. |0J). This value is 
employed in many works (see, e.g. Refs. 0, 0, Q, 0, 
01). The rjL=i values used in calculations of the main 
properties of the ISGDR in nuclei under consideration 
are given in Tabled 

The isoscalar monopole and quadrupole GRs in nuclei 
under consideration are experimentally studied in many 
works Q S H m El 113 . The deduced strength distri- 
butions are presented either in terms of the peak energy 
^peak total width (obtained by the Lorentzian 
or Gaussian fit), or in terms of the mean energy uj^ and 
RMS energy dispersion (obtained for a certain exci- 
tation energy interval). These data are used below for 
comparison with the results, obtained within the semi- 
microscpopic approach. The energy-averaged strength 
fimctions S'l(w) calculated with use of the radial part 
^Lir) = {L — 0,2) of the first-order probing op- 
erators exhibit well prominent resonances in their en- 
ergy dependence. The respective peak energies and to- 
tal widths are given in Table together with the cor- 
responding experimental values. The relative energy- 
weighted strength functions ?/l(w) = ujSl{uj)/ {EWSR)l 
are shown in Fig. la (L = 0) and in Fig. 2a {L = 2) 



FIG. 3d: The reduced energy-dependent transition density 
calculated at the peak energy of the LE-ISGDR (Table^} . The 
notations are the same as in FIG. la. Curves with filled circles 
and open circles correspond, respectively, to the additional 
maximuma for ^*Ni and ■^°*Pb as given in Table Q 




m (MeV) 



FIG. 4: Differential partial proton branching ratios for the 
HE-ISGDR in ^'^*Pb calculated for several decay channels 
(Sfj, — 1). The thick, thin, dashed, dotted and dash-dotted 
lines are for direct decay into 3si/2, '2d^/2, 1^*11/2, 2d5/2, and 
I37/2 one-hole states of ^"^Tl, respectively. 



for nuclei under consideration. Microscopically the IS- 
GMR and ISGQR are mainly due to 2fujj particlc-hole- 
type excitations. The low-energy component of the IS- 
GQR (Fig. 2a) is due to Qhto excitations corresponding 
to single-particle transitions with changing both the ra- 
dial and orbital quantum numbers and having, therefore, 
a small relative strength. The calculated strength func- 
tions are used to evaluate (for a certain energy interval) 
the parameters of the ISGMR and ISGQR: mean energy 
ujl, RMS energy dispersion A^, relative strength xl- 
These values are listed in Table Ull Comparison with the 
available experimental data is given in Tables Hvl and Ivl 
The reduced transition densities RL{r,ojp^^f.) calculated 
for the main-tone isoscalar monopole and quadrupole 
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FIG. 5a: The calculated relative energy- weighted strength 
function for the ISGMR2. The notations are the same as in 
FIG. la. 



FIG. 5b: The calculated relative energy-weighted strength 
function for the ISGQR2. The notations are the same as in 
FIG. la. 



GRs in nuclei under consideration are shown in Fig. lb 
and Fig. 2b, respectively. The transition densities are 
used to evaluate parameters ?]l accordingly to Eq. 
These values are given in Tabled 



B. Properties of the ISGDR 

The gross properties of the ISGDR in nuclei under 
consideration are described within the approach using 
the isoscalar second-order dipole probing operator with 
parameters ?7l=i taken from Table H The calculated 

~ (2) 

energy-averaged strength function S\^J^^{uj) (Fig. 3a) ex- 
hibits a "bi-modal" energy dependence, corresponding 
to the low- and high-energy components of the ISGDR 
(LE- and HE-ISGDR, respectively). As a resonance in 
the strength function energy dependence, the HE-ISGDR 
is well prominent and can be described in terms of the 
peak energy and total width, while the LE-ISGDR is 
less prominent. The peak energies of both components 
are given in Table |3 together with the latest experimen- 
tal data. The relative energy-weighted strength func- 
tions ?yr 1,1 (iv) are very close to those, shown in Fig. 1 of 
Ref. [i3 for nuclei under consideration (except for ^^Ni) 
and are not shown here. The parameters of both compo- 
nents calculated with the use of the above strength func- 
tions are given in Table lllll A part of these results related 
to certain excitation-energy intervals are compared with 
the available experimental data deduced for the same in- 
tervals f Table HVjl . The strength function yfli{i-^) cal- 
culated for ^*Ni is shown in Fig. 3b. Because the main- 
tone transition density is nodeless, the energy-dependent 

(2) 

transition density R)^!^^{r,ujpeak) calculated at the peak 
energy of each component exhibits one-node radial de- 
pendence (Figs. 3c, 3d). 

Direct nucleon decay of giant resonances is closely re- 
lated to their microscopic structure. For this reason, the 



TABLE VI: Calculated partial branching ratios for direct 
proton decay of the HE-ISGRs in ^**Ni (S^ = 1). Branching 
ratios are given in %. 
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35.0 
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decay probabilities belong to the main properties of GRs 
together with the energy, total width, and transition den- 
sity. The use of the CRPA method together with the phe- 
nomenological treatment of the spreading effect allows 
us to evaluate the partial direct-nucleon-decay branching 
ratios for various GRs within the semi-microscopical ap- 
proach. Turning to the HE-ISGDR in nuclei under con- 
sideration, we present the direct-nucleon-decay branch- 
ing ratios calculated with the use of unit spectro- 
scopic factors Sf^ for single-hole states in closed-shell sub- 
systems and of unit ratio S^/vl for single-quasiparticle 
states in open-shell subsystems iTables lVllXII . Some par- 
tial branching ratios corresponding to population of deep- 
hole states are not shown, but they are included in the 
values of the respective total branching ratios also given 
in the tables. The recent experimental data of Refs. j^Q 
on the partial direct-proton-decay branching ratios for 
the HE-ISGDR in ^ospb are given in Table EH together 
with the respective calculated values obtained with the 
use of the experimental spectroscopic factors. Calculated 
for the same resonance, the differential partial proton 
branching ratios of Eq. (|3J) are shown in Fig. 4. One can 
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FIG. 6a: The reduced energy-dependent transition density 
calculated at the peak energy of the ISGMR2 (Table |lj. The 
notations are the same as in FIG. la. 
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FIG. 7: The calculated relative strength function for isoscalar 
GRs in ^°*Pb. The thick, thin, dashed, dotted, and dash- 
dotted lines are for the ISGMR, ISGQR, ISGDR, ISGMR2, 
and ISGQR2, respectively. 
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FIG. 6b: The reduced energy-dependent transition density 
calculated at the peak energy of the ISGQR2 (Table QJ. The 
notations are the same as in FIG. la. 

see from this figure the role of the penetrability factor 
in formation of the HE-ISGDR in different proton decay 
channels. 



C. Main properties of the ISGMR2 and ISGQR2 

We describe the gross properties of the isoscalar 
monopole and quadrupole overtones in nuclei under con- 
sideration, using the respective second-order probing op- 
erators with parameters tjl {L = 0, 2) taken from Ta- 
ble n The overtone strength is distributed over a wide 
energy interval exhibiting the main peak at a high ex- 
citation energy. The w^g^j. values are given in Table 

The relative energy-weighted strength functions y)^ (w) 
are shown in Figs. 5a (L = 0) and 5b (L = 2), while 
the parameters of the high-energy components are given 



in Table Ixlll Because the main-tone transition density 
has one node inside of the nucleus (L = 0, Fig. la) or is 
nodeless {L = 2, Fig. lb), the overtone transition den- 
sity has two nodes {L = 0, Fig. 6a) or one node (L = 2, 
Fig. 6b), respectively. The direct-nucleon-decay branch- 
ing ratios calculated for the high-energy components 
of both overtones with the use of unit spectroscopic fac- 
tors, are given in Tables IVIIXI Partial branching ratios 
for direct proton decay of the above resonances in ^ospj-, 
are also calculated with the use of experimental spectro- 
scopic factors (Table Drr|) . 



IV. DISCUSSION OF RESULTS AND 
SUMMARY 

Within the CRPA-based semi-microscopic approach, 
all the main properties of a given GR can be described 
in a transparent and rather simple way using universal 
parameters for medium-heavy and heavy (spherical) nu- 
clei. The symmetries of the model Hamiltonian are re- 
stored via the respective selfconsistency conditions. In 
the present version of the approach the spin-orbit part of 
the nuclear mean field is mainly responsible for incom- 
plete restoration of translation invariance of the model. 
If this part was taken equal to zero, the SS would ex- 
haust more than 99.5% of the respective sum rule, while 
the gross properties of isoscalar GRs are practically not 
changed. The selfconsistency of the present version of 
the approach can apparently be improved, provided that 
the isoscalar spin-dependent part of the Landau-Migdal 
particle- hole interaction is taken into account. The re- 
spective study is outside the scope of the present work 
and will be addressed in a future publication. 

Discussing the results, we start from the gross proper- 
ties of the isoscalar GRs. Taking ^°^Pb as an example, 
one can see from the results of the semi-microscopic cal- 
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TABLE VII: Calculated partial branching ratios for direct 
nucleon decay of the HE-ISGRs in ^"Zr (S,! — 1). Branching 
ratios are given in %. 
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culations of the relative energy-weighted strength func- 
tion yhi'^) (Fig. 7) the general tendency for changing 
isoscalar strength distribution with increasing excitation 
energy. The main components of the ISGMR2 and IS- 
GQR2 arc not well-collectivized exhausting a not-too- 
large part of the respective total strength (Table IXII|I . 
As a result, the A-dependence of cj^f^^. for these over- 
tones is not so regular, as it takes place for more col- 
lective GRs (Table . According to the data shown in 
Tables HlllVl 13 the experimental energies of the ISGMR, 
ISGQR, and HE-ISGDR in nuclei from a wide mass in- 
terval are reasonably described within the approach. The 
energy of the LE-ISGDR in the same nuclei is described 
satisfactorily (Table HIIIVl). The energy of the recently 
found ISGQR2 in ^ospb 26.9 ± 0.7 MeV 0, is also 
satisfactorily reproduced in the calculations (Table QJ. 
Thus, the calculated energies of the second isoscalar gi- 
ant quadrupole and monopole resonances in nuclei under 
consideration can be used as a guide in search of these 
GRs experimentally. 

We try to elucidate the universal phenomenological de- 
scription of the total width of an arbitrary GR using an 
appropriate smearing parameter with the saturation-like 
energy dependence of Eq. ||2Jl. In applying to isovec- 
tor GRs, a similar attempt was found to be satisfac- 
tory "2^. According to this description, the total width 
of low-energy GRs (except for the isobaric analog res- 
onance) is mainly due to the spreading effect. Such a 
case is realized for the main-tone isoscalar monopole and 
quadrupole resonances. As follows from CRPA calcula- 
tions, a significant part of the total width of high-energy 
GRs is due to the particle-hole strength distribution and 



TABLE VIII: Calculated partial branching ratios for direct 
nucleon decay of the HE-ISGRs in "''Sn {S^ = 1 and S^, = 
are used in calculations of proton and neutron branching 
ratios, respectively). Branching ratios are given in %. 
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coupling to the continuum. The rest is due to the spread- 
ing effect, which leads also to averaging the strength dis- 
tribution over the energy. Regarding the experimental 
total widths of the isoscalar GRs in nuclei under con- 
sideration, we note that the widths of the ISGMR and 
ISGQR, the RMS energy dispersion for both components 
of the ISGDR are well described within the approach (Ta- 
bles |l||V|). The experimental total width of the ISGQR2 
in 208pb 6.0 ± 1.3 MeV 5] agrees well with the RMS 
energy dispersion, calculated for the main component of 
this resonance (Table lXIlj) . It is worth to note the notice- 
able scatter in the experimental data on the total width 
of the HE-ISGDR 3, 5, 6]. This scatter and also the scat- 
ter of experimental relative strengths xl are apparently 
explained by uncertainties in the substraction of the un- 
derlying continuum in the analysis of the (a, a')-reaction 
cross section. For this reason we do not compare in the 
present work the calculated xl values (Tables Hll and ITTH) 
with the respective experimental data. Nevertheless, one 
can compare the calculated (Fig. 3b) and experimental 

(Ref. 0) strength functions y^lii^^) for ^*Ni and dis- 
cover satisfactory agreement. 

As expected, the radial dependence of the one-node 
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TABLE IX: Calculated partial branching ratios for direct 
neutron decay of the HE-ISGRs in ^''^Sm [S^ = 1). Branching 
ratios are given in %. 





(15-35 MeV) 


(21-38 MeV) 


(24-36 MeV) 


(1/2) + 


1.9 


2.0 


1.0 


(3/2)+ 


3.4 


3.3 


2.2 


(11/2)- 


5.2 


4.9 


9.9 


(5/2)+ 


5.1 


5.7 


4.4 


(7/2)+ 


2.6 


4.1 


3.7 


(9/2) + 


1.8 


4.3 


5.2 




23.1 


35.1 


35.3 



transition density RL{r,ujpeak) for the ISGMR and HE- 
ISGDR (Figs, lb and 3c) is rather close to that of the 
respective transition density calculated within the scaling 
model 32], provided the ground-state density of Eq. Q is 
used in the calculations. It can be seen, for instance, from 
Fig. 2 of Ref. 18]. We note also, that the difference of the 
transition densities calculated at the peak energy of each 
ISGDR component is not so large (Figs. 3b and 3c) to 
allow for searching for exotic explanations of LE-ISGDR 
nature. As it is also expected, the overtone transition 
density has one extra node inside the nucleus relative to 
the main-tone transition density (Figs, lb, 2b, 3c, 5b, 6b). 
Concluding consideration of the gross properties of the 
isoscalar GRs, we note that the use of the modified 
smearing procedure (as compared with that of Ref-fTsf) 
leads practically to the same results. It can be seen from 
comparison of the data from Tables ITU and IIIl| with those 
from Tables I and II of Ref. ; respectively. 

The use of the modified smearing procedure (Sub- 
sect. Ill B|) allows us to describe within the present 
approach the direct-nucleon-decay branching ratios for 
high-energy GRs. The respective calculation results ob- 
tained for the isoscalar overtones in nuclei under con- 
sideration, are given in Tables IVIIXI For the ISGDR, 
the total direct-nucleon-decay branching ratio decreases 
with decrease of the peak energy (and with increase of the 
mass number) . The relative change is larger for the total 
proton branching ratios due to the difference in the pen- 
etrability factors. We note, that the results of Ref. |19l ]. 
where the main shortcoming in evaluation of the direct- 
nucleon-decay branching ratios for high-energy GRs was 
eliminated (Subsect. lII B|) . are close to those of Tables lVTl 

nil 

The recent experimental data of Refs. 0,0] on partial 
direct-proton-decay branching ratios for the HE-ISGDR 
in •^"^Pb are satisfactorily described within the present 
approach, provided that experimental spectroscopic fac- 
tors for the final single-hole states of 2°'^T1 (Table ^ 
are taken into account. In Ref. direct neutron de- 
cay of the same resonance into the final states of ^°^Pb 
from excitation-energy interval 0-6 MeV has been also 
observed. The deduced branching ratio 23 ± 5% reason- 



TABLE X: Calculated partial branching ratios for direct nu- 
cleon decay of the HE-IS CRs in ^"^^Pb (S^ = 1). Branching 
ratios are given in %. 
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TABLE XI: Partial branching ratios for direct proton decay 
of the HE-ISGRs in ^°^Pb into some one-hole states of ^°'^T1. 
Experimental spectroscopic factors taken from Ref. |3lll 
are used in calculations. Excitation-energy intervals are taken 
the same as in Table ^ Calculation results for decays of 
the HE-ISGDR are compared with the experimental data of 
Refs. 6, 7]. Branching ratios are given in %. 
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'Preliminary results 



ably agrees with the value 15.2% obtained with the use 
of unit spectroscopic factor for the respective one-hole 
states (Table 1^. Reasonable description of the above 
data allows us to infer that the calculation results shown 
in Tables IVIIIXI will be useful for the analysis of forth- 
coming experimental data on direct nucleon decays of 
the HE-ISGDR in several medium-heavy and heavy nu- 
clei 8]. Some evidence for direct proton decay of the IS- 
GQR2 in ^°*Pb in the coincidence (a, a'p) experiments 
have been recently found 0, 0] . It allows us to hope that 
the calculated parameters of the ISGQR2 and ISGMR2 
(including the energy, total width, transition density, and 
branching ratios for direct-nucleon-decay) would be also 
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TABLE XIL Paramaters of the ISGMR2 and ISGQR2 cal- 
culated for a certain excitation-energy interval. All the pa- 
rameters are given in MeV except for x, which is given in 
%. 
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useful for experimental search of these resonances. 

In conclusion, we extend a CRPA-based partially self- 
consistent semi-microscopic approach to describe direct 
nucleon decay of high-energy giant resonances. The main 
properties of the isoscalar overtones (ISGDR, ISGMR2, 
ISGQR2) in a few singly- and doubly-closed-shell nuclei 
are described within the approach and found to be in 



reasonable agreement with available experimental data 
including the latest ones. Abilities of the approach for 
description of the main-tone resonances are successfully 
checked. Predictions concerning forthcoming experimen- 
tal data on the isoscalar overtones are also presented. 
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